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ABSTRACT 

Galaxy growth depends critically on the interplay between radiative cooling of cosmic gas and the 
resulting energetic feedback that cooling triggers. This interplay has proven exceedingly difficult 
to model, even with large supercomputer simulations, because of its complexity. Nevertheless, real 
galaxies are observed to obey simple scaling relations among their primary observable characteristics. 
Here we show that a generic emergent property of the interplay between cooling and feedback can 
explain the observed scaling relationships between a galaxy’s stellar mass, its total mass, and its 
chemical enrichment level, as well as the relationship between the average orbital velocity of its 
stars and the mass of its central black hole. These relationships naturally result from any feedback 
mechanism that strongly heats a galaxy’s circumgalactic gas in response to precipitation of colder 
clouds out of that gas, because feedback then suspends the gas in a marginally precipitating state. 
Subject headings: galaxies: evolution — galaxies: star formation 


1. INTRODUCTION 

Radiative cooling is essential to galaxy formation be¬ 
cause shock waves heat cosmic gas to million-degree tem¬ 
peratures as gravity pulls it into galaxies. Early mod¬ 
els of galaxy formation therefore focused on comparing 
the cooling time tcooi of the gas with the free-fall time 


ation Q (Rees & Ostriker 

1977 Silk 

1977 

Binney 

1977 

& Kees 

1985|). Gas with Gooi < tff was presumed 


of a multiphase medium consisting of cold star-forming 
clouds embedded in an ambient million-degree atmo¬ 
sphere. Later models recognized that uncompensated 
cooling would produce many more stars than are ob- 
serv ed, requiring a feedback response that inh i bits cool- 


i:^’renk|l991 

in small ga 

). Supernovae can provide suthcient feedback 
laxies but cannot fully compensate for cool- 

ing in larger ones 

Benson et al. 12003 

Groton et al. 12006 

Bower et al.|2006 

. Instead, accretion of cooling gas onto 

a central supermassive black hole seems required to gen¬ 
erate the feedback that regulates star formation in the 

universe’s 1 

argest galaxies (Tabor & Binnev|1993 Silk & 


Exactly how black-hole feedback becomes tuned to 
compensate for radiative cooling at distances > 10^ par¬ 
secs has been a resilient mystery, but there have been 
some recent breakthroughs. It now appears that the 
black-hole fueling rate in giant galaxies depends sensi- 


tively on development of a multiphase medium (Pizzo- 
lato fc Sokerj 2005 |Voit et ahl 2008[ |Cavagnolo et al. 


2008p. Numerical sirnulations indicate that circumgaiac- 


tic gas in approximate thermal balance develops a mul- 


cipitate out of the hot medium when Gooi drops below 
~ lOtff and then accrete chaotically onto the galaxy’s 


(jGaspari et al. 

2012 

2013 

2014 

Li & Bryan 

2014a|b 


galactic medium, raising tcooi and diminishing the pre¬ 
cipitation. Interplay between precipitation and feedback 
therefore causes min(tcooi/^ff) to fluctuate in the range 
5 < min(tcooi/tff) ^ 20 and to settle near Aooi/tff ~ 10. 

X-ray observations of circumgalactic gas around gi¬ 
ant galaxies are revealing a distinct cooling-time floor a t 
tcooi ~ lOtff ( Voit fc Donahue|[2015 Voit et aL]|20 
confirming a key prediction of those numerical mode 


b) 


(Sharma et al. 2012aI. Here we explore the implications 


(IPizzolato & Soker 

2010 

McGourt et al. 

2012 

Sharma 

, , 3kT 

hne +ni'\ 
^ 2ni y 

et al. 

2L)12bl. in simulated galaxies, co 

d clouds pre- 

lOAi, 
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lor galaxy evolution if feedback driven by precipitation 
of cold clouds produces a similar cooling-time floor in 
galactic systems of all sizes, from large galaxy clusters 
all the way down to the smallest dwarf galaxies, and we 
show that this simple hypothesis can account for many 
of the observed properties of galaxies. 

2. THE PRECIPITATION THRESHOLD 

There is not yet a complete theory for the origin of the 
critical Gooi/tfr value. In lieu of a complete theory, we 
follow the observations and assume that Aooi/tff ~ 10 
is the threshold at which precipitation triggers a strong 
feedback response that reheats the circumgalactic gas. 
If the amount of precipitating gas rapidly increases as 
tcooi/tff drops, then any form of feedback fueled by pre¬ 
cipitation causes the system to linger near Gooi/te ~ 10, 
regardless of whether the feedback comes from black-hole 
accretion or supernova explosions. The electron number 
density at radius r in a generic precipitation-regulated 
system is therefore 


( 1 ) 


where T is gas temperature, Ui i s ion number density , 
and A is the usual cooling funtion ( Sharma et al.|20i^ ). 
In such a system, the radius r^oo that encompasses an 
average mass density 500 times the critical density is an 
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approximate upper limit on the size of the precipitating 
region, because fcooi ~ at that radius. Intriguingly, 
spectroscopic absorption-line studies of the circumgalac- 
tic medium indicate that small clouds of low-ionization 
gas are com mon inside of rsnn bu t are rare outside of 
that radius (Liang & Chen 2014), suggesting that the 
low-ionizatioh material might be incipient precipitation. 

Condensation and precipitation of the circumgalactic 
medium supply the cold gas needed for star formation 
in a precipitation-regulated galaxy. One expects that 
gas supply to be proportional to the integral of 
within the volume bounded by rsoo. We therefore express 
the time-averaged star formation rate in a precipitation- 
regulated galaxy as 


M* 


Stt GfigmpkT 

M M 


e*/b-^500 


( 2 ) 


where e* is an efficiency parameter accounting for the 
complexities of star formation and feedback, /b is the 
cosmic baryon mass fraction, /ierup is the mean mass 
per electron, M 500 is the total mass within rsoo, and the 
approximation assumes that T and r/tg are independent 
of radius. (In giant galaxies, the limiting value of M* 
is smaller because the cosmological upper limit on Ue 
is more stringent at large radii than the precipitation 
threshold.) 

According to equation (1), feedback in a precipitation- 
regulated galaxy limits cooling by driving down the den¬ 
sity of its circumgalactic medium. It does not perma¬ 
nently expel all the circumgalactic gas from the system 
but rather inflates it like a balloon, pushing gas that 
would have been inside of r^oo to larger radii, thereby 
lowering the mean density of the precipitating region. 
Electron densities around dwarf galaxies, which have 
comparatively low T and large A, are consequently sev¬ 
eral orders of magnitude smaller than those around giant 
galaxies. 

According to equation (2), the average star-formation 
rate of a precipitation-regulated galaxy depends inversely 
on chemical enrichment of the circumgalactic gas through 
the cooling function A, which is proportional to the abun¬ 
dances of elements heavier than H and He for gas in 
the 10® — 10^ K temperature range. This inverse depen¬ 
dence arises because enrichment makes the circumgalac¬ 
tic medium more volatile. Paradoxically, adding heavy 
elements to a precipitating system reduces M* because 
feedback is triggered at a lower circumgalactic gas den¬ 
sity. Observations of real star-forming galaxies show that 
they share this qualitative feature, in that their star- 
formation rates anticorrelate with chemic al enrichment 
( Ellison et al.|[2008 Mannucci et al.|[2010 |. 


3. SATURATION OF ENRICHMENT 


Coupling between chemical enrichment and star for¬ 
mation in a precipitation-regulated galaxy causes its en¬ 
richment level to saturate. In order to illustrate how 
saturation arises, we adopt a simplistic assumption: that 
the heavy elements made by stars mix evenly with all the 
gas associated with a galaxy, including its circumgalactic 
medium. The mass Mz,gas of heavy elements in the gas 
then accumulates according to Mz,gas = {Y — ^gas)Af*, 
where Y is the heavy-element yield from a single gener¬ 
ation of star formation and Z^as = Mz,gas/Mgas- Note 


that Mgas = fhM^oo — M* represents the entire gas mass 
originally associated with the galaxy’s dark matter, less 
the amount that has turned into stars, and generally ex¬ 
ceeds the gas mass within 

This chemical enrichment model is oversimplified but 
nicely clarifies the mechanism of saturation. Taking the 
derivative of Zgas with respect to time, we find 


Zgas — 


EM* 

ZgasfhM^OO 


1 

Mgas 


( 3 ) 


The first term in brackets is a source term for element 
production, and the second is a dilution term account¬ 
ing for inflowing unenriched gas. Together they tune 
both the galaxy’s heavy-element abundance and its star- 
formation rate to particular values determined by the 
galaxy’s gravitational potential well. Early in a galaxy’s 
history, when enrichment is low, the source term domi¬ 
nates and Zgas increases. However, enrichment also in¬ 
creases the ability of circumgalactic gas to cool, which 
promotes precipitation. Feedback therefore causes the 
circumgalactic gas to expand until it reaches the pre¬ 
cipitation threshold at tcooi/te « 10. In this state of 
marginal stability, expansion continues while Zgas in¬ 
creases, which drives down the star-format ion rate. As a 
result, the heavy-element abundance of a precipitation- 
regulated galaxy plateaus when it reaches saturation at 
Z«EM*(/bM5oo)-i. 

Chemical enrichment of real galaxies is undoubtedly 
more complicated. However, heavy-element abundances 
in any precipitation-regulated system will still saturate 
as long as enrichment of the circumgalactic medium is 
proportional to the total production of elements by the 
galaxy’s stars. In such a system, chemical enrichment of 
the circumgalactic medium reduces its condensation rate 
at the precipitation threshold, which is what is required 
for saturation to happen. 


4. GALACTIC STAR-FORMATION HISTORIES 

Approximate star-format ion histories for precipitation- 
r^ulated galaxies can be obtained by rewriting equation 
dW as 


H 


^^Ml 
dt * 


(2 X 10 “ 


Y 


riQ i;2oo 


( 4 ) 


where H is the Hubble expansion parameter and ti 2 oo 
is the circular velocity Vc ~ in units of 

200kms“^. To get equation (|3, we assumed A « 
(1.8 X 10~^^ ergcm ® s~^)(Z/Zf 7 ))(T/ 10 ® K)~\ wh ich is 
within 50% of the [Sutherland fc Dopita (1993 1 tab¬ 
ulation for 10® K < T < 10 ' K. We also assumed 
Z « EM*/(/bM 5 oo) and kT « ^iripV^, as appropriate 
for hydrostatic gas with Ue oc r~^ in an isothermal po¬ 
tential. Mass growth of a typical galactic halo follows 
M 500 = with 0 ss 0.5(gfit)(l 


zl^®® m a 


ACDM cosmology (Neistein & Dekel 200_8|, and using 
this approximation to integrate equation gives 

1/2 

(5) 


M* « (2 X 10'^^Mq) 


^0 

Y 


where 


F(t,t{) = 


Hi 


,10 


it') 




dt' 


( 6 ) 
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Fig. 1.— Specific star-formation rates of galaxies over cosmic 
time. Points with error bars show observations of specific star- 
for mation rates in galaxie s compiled from many different sources 
by | |Behroozi et al.|2013b[ |, color-coded according to galactic stel¬ 
lar mass: < 10 '='Mq (purple), (magenta), 

(brown), (gold), (orange), > IOIIM0 

(red). Each point represents an average over many galaxies in that 
stellar-mass bin at a particular epoch in cosmic time, and the error 
bars indicate uncertainty in the mean. The green solid line shows 
specific star-formation rate predictions for t{ = 0. Dashed lines 
show predictions for if = 2, 4, 6, and 8 Gyr. The dotted line shows 


and t{ is a “formation time” at which star formation 
begins. 

Star-formation histories derived from equation ([^ 
align with measurements of specific star-format ion rates 
in galaxies and how they trend across cosmic time (Fig¬ 
ure!^. In this context, the precipitation thr eshold be¬ 
haves like a regulator of the sort envisioned by [Lilly et al.| 
(2013|). The circumgalactic medium acts as a reservoir 
that collects both low-enrichment gas accreting onto the 
galaxy and highly enriched gas flowing out of the galaxy. 
Precipitation is how gas goes from the reservoir into the 
galaxy, where it can form stars, and the gas transfer rate 
depends on T and A. The circumgalactic gas temper¬ 
ature depends in turn on the depth of the galaxy’s po¬ 
tential well, while A depends also on enrichment of the 
circumgalactic medium, which self-regulates through the 
mechanism of abundance saturation. 


5. ENRICHMENT AND STELLAR MASS 

Precipitation-regulated galaxies exhibit scaling rela¬ 
tions among galactic stellar mass, enrichment level, and 
total mass that are nearly identical to the observed rela¬ 
tions. The mean enrichment level is Z « FAL*/(/bM 5 oo), 
which corresponds to 


Z 


0.6 



0.7 

(e*F)0-3 


IQii Mq 


0.4 


w ’ 

tlQ 


according to equation (5). The normalization of this re¬ 
lation is similar to observations but depends somewhat 
on the uncertain parameter e*, which requires more de¬ 
tailed modeling. Instead, the scaling with M, is what 
provides the most compelling evidence in support of the 
precipitation-regulated framework. 

Figure shows how that scaling corresponds to ob¬ 
served ZpVf*) relations. Precipitation-driven feedback 


Fig. 2.— Relationship between total stellar mass M* and heavy- 
element abundance Z. The maroon line shows how Z{Mt,) scales 
for precipitation-regulated galaxies; the grey line shows the approx¬ 
imation Z ^ (ALt/lO^^M©)'^-^. Bl ue points show stellm Fe abun¬ 
dances in low-mass galaxies from [McConnachiej (|2012|l. Orange- 
and salmon-colored regions show stella r Mg and he abunda nces, 
respectively, in elliptical galaxies from [Graves et al.| ||2009||. In 
green is the mean gas-phase abund ance derived by Kewley fe El-j 
jlisonj (|2008|l from data compiled by [Tremonti et al. (|2UU4||, along 
with Iff and 2cr standard-deviation contours. 


naturally leads to heavy-elem ent abundances in dwarf 
galaxies (McConnachie 2012) tha t are only ~ 1% o f 
the levels found in large galaxies ([Graves et al.||2009|), 
with little dependence on environnient. INIo dependence 
of the parameter e* on galaxy mass is required. Ra¬ 
diative cooling is simply more efficient at 10® K than 
at 10^ K, making the gas around smaller galaxies more 
susceptible to precipitation. As a result, their circum¬ 
galactic gas densities are lower, and star formation lim¬ 
ited by the flow of condensing clouds into the galaxy is 
correspondingly slower, which strongly limits the rate of 
heavy-element enrichment. Furthermore, the model pre¬ 
dicts Z oc for late-blooming galaxies 

lying above the green line in Figure in quantitative 
agreement with observatio ns of galaxies with large spe¬ 
cific star-formation rates (Mannucci et al. 2010 Salim 


et al. 


2014), since M*/M* cc F ® for such galaxies. 


This regulation mechanism also accounts for the ob¬ 
served scaling of stellar baryon fraction with Vc (see Fig¬ 
ure , predicting that 


f* = 




fhM^OO 




Y 


Hn 


^200 > 


( 8 ) 


for present-day galaxies, for which F{t,0) « 0.44. An¬ 
other way to characterize this dependence is in terms of 
the well-k nown Faber-Jackson and Tully-Fisher scalin g 
relations ([Faber & Jackson||1976| [Tully & Fisher||1977|), 
which show that the stellar luminosity L of a galaxy de- 
pends on its gravitational potential according to L cx 
with ^ « 4 — 5. A precipitation-regulated galaxy should 
have M* cx u®, which gives L oc u® for an invariant stel¬ 
lar mass-to-light ratio. Recent observations show that 
the stellar mass-to-light ratio for massi ve galaxies de¬ 
pend s on galactic mass, with M*/L oc Vc (Spiniello et al. 
2014), which leads to L oc for massive galaxies and 
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Fig. 3. — Relationship between a galaxy’s stellar baryon fraction 
/* and the maximum circular velocity Uc,max of its gravitational 
potential. The maroon line shows the dependence of /* on Uc,max 
for precipitation-regulated galaxies; the grey line shows the approx¬ 
imation /* a:i 0.25u^nn- Points show measurements of /* in spiral 
galaxies from |McGaugh| ( |2005^ . 

indicates that the Faber-Jackson and Tully-Fisher rela¬ 
tions are natural outcomes of precipitation-regulated star 
formation. 

6. GIANT GALAXIES 


Fig. 4.— Relationship between velocity dispersion cr„ of a 
galaxy’s stars and the mass Mbh of its central black hole. Data 
points for ellipticals (circles) , classical bulges (squa res), and pseu¬ 
dobulges (triangles) are from|Kormendy & Ho|||2013||. Star symbols 
show how galaxies belonging to galaxy clusters sliift to the right 

when an effective velocity dispersion (t„ gfl ~ 

is used. Dashed green lines show model predictions for different 
values of the efficiency factor cbu- A solid green line shows the 
lower limit on ALbh derived for cbh = 0-1. 

The cooling time at this boundary in massive galaxy clus¬ 
ters is 


The trends shown in Figures and must change in 
giant galaxies for two reasons: First, giant precipitation- 
regulated galaxies with deep potential wells transform 
their gas into stars more quickly than smaller galaxies, 
producing the trend in Figure which must stop before 
/* reaches unity. According to equation (|^, galaxies 
with the greatest values of Vc are the first to become 
gas-limited. This event inevitably causes star formation 
to decline at a cosmic time approximately proportional 
to v~'^, possibly explaining the phenomenon astronomers 
call downsizing (Madau & Dickinson 2014). Second, as 
the baryon balahce tips from gas toward stars, black- 
hole feedback and supernova explosions produced by a 
giant galaxy’s old stellar population become capable of 
shutting d own precipitation and preventing further star 
formation (|Voit et al.|2015b|). Observations in dicate that 
galaxies reach this threshold when /* ~ 0.3 (McGaugh 


et al. 2010 Dai_et al. 2010). According to the trend 
shown in Figure [s] galaxies with Vc ~ 250kms“^ are the 
ones currently making this transition. The corresponding 
stellar mass at the transition is ~ 10^^Mq and should 
change with time oc in accordance with the observed 
lack of time dependenc e in this transitional stellar mass 
(|Behroozi et al.||2013a|). 

Further increases m total mass cause a qualitative 
change in how cosmic systems acquire stars. Once the 
system temperature exceeds ~ 1 keV, the precipita¬ 
tion model predicts that /* should decline with increas¬ 
ing system mass, because the precipitating region starts 
to shrink relative to ^ 500 . Its outer boundary in this 
regime is given by the radius rprecip at which the elec¬ 
tron density at the precipitation threshold equals the 
value fbi^irG, obtained by assuming that the 
electron density is proportional to the total mass density. 


t f(r (9) 

tcooUrprecipj ~ lyj 

Myr. 

Consequently, the mass cooling rate within the precipi¬ 
tation zone scales as M* oc e*f;^(rprecip), implying that 
the stellar mass of the central galaxy decouples from the 
ambient gas temperature and stops growing in propor¬ 
tion to the total mass of the system. High-mass systems 
such as galaxy clusters must therefore acquire the ma¬ 
jority of their stars by cannibalizing lower-mass systems 
rather than by producing new stars themselves. 


7. CENTRAL BLACK-HOLE MASS 

These mass-dependent features of the precipitation- 
regulated mo del appear to be reflect ed in the observed 
relationship (Kormendy & Ho 2013) between the mass 
Mbh of a galaxy's central black hole and the line-of- 
sight velocity dispersion cry « Vc of its stars. If mass 
accretion onto a central black hole is the primary energy 
source for regulating precipitation, then a galaxy’s black- 
hole mass Mbh should depend on the cooling rate of its 
precipitation zone via 


cbhI^c^Mbh 


J Airr^nlAdr 


( 10 ) 


3/b 

bG^errip 


kTal min 


^ ^ff(F50o) 

^ff(^precip) 


where cbh is a coupling efficiency representing the frac¬ 
tion of accreted mass-energy that ends up as heat in 
the circumgalactic medium. In massive galaxies with 
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cr„ ^ 250kms“^, the model predicts Mbh oc ct® and 
agrees with the data for esH ~ 10“^ — 10“^ (Figure 0. 
In precipitation-regulated galaxies with (7„ < 200 km s^ , 
the model predicts a steeper relation (Mbh oc crj,), which 
agrees with the data for a similar coupling efficiency. At 
the high-mass end, the model predicts Mbh cx Ta^, 
which explains why the masses of black holes in the 
central galaxies of galaxy clusters tend to exceed those 
in other large galaxies with similar velocity dispersions. 
One can account for the decoupling of central cluster 
galaxies from system temperature by defining an effec¬ 
tive velocity dispersion cr^ oc Tcrj^. Plotting Mbh 
as a function of cr^^eff instead of cr„ then shifts central 
cluster galaxies onto a tighter relation corresponding to 
cbh « 10-2-3. 


8. TOWARD GREATER REALISM 

Real galaxies are clearly more complicated than por¬ 
trayed in these spartan precipitation-regulated models. 
For example, the models completely ignore angular mo¬ 
mentum, which is critical for the formation of galactic 
disks and is likely to shi ft the precipitation t hreshold to 
lower values of tcooi/iff (|Gaspari et al. 2014). They do 
not account for cold clouds that tall into galaxies with¬ 
out being shock-heated to the ambient gas temperature. 
They do not attempt to describe the enrichment gradi¬ 
ents observed in real galaxies. Nor do they attempt to 
account for environmental effects, such as stripping of 
the circumgalactic medium from galaxies orbiting within 
larger gravitationally-bound systems. 

On the other hand, the models do a remarkable job of 
reproducing several of the most important scaling rela¬ 


tions observed among global properties of galaxies, par¬ 
ticularly since they rely on just a few highly plausible 
assumptions: (1) that galaxies acquire star-forming gas 
through condensation of the circumgalactic medium, (2) 
that condensation triggers feedback which heats the cir¬ 
cumgalactic medium, and (3) that chemical enrichment 
of circumgalactic gas is proportional to the amount of 
star formation in the galaxy. The precipitation thresh¬ 
old therefore seems likely to be an important feature of 
galaxy evolution, not just in the giant galaxies where 
it is directly observed but also in lower-mass galaxies, 
around which the volume-filling gas component is essen¬ 
tially unobservable. Such a threshold should be an emer¬ 
gent feature of any numerical simulation of galaxy evolu¬ 
tion in which condensation of circumgalactic gas triggers 
strong feedback that suspends the ambient medium in a 
marginally precipitating state. If that happens, then star 
formation fueled by precipitating gas will lead to satura¬ 
tion of chemical enrichment and will produce populations 
of galaxies with mass-dependent trends strikingly similar 
to those observed among real galaxies. 
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